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Abstract 

We consider a generic optomechanical system, consisting of a driven 
optical cavity and a movable mirror attached to a cantilever. Systems of 
this kind (and analogues) have been realized in many recent experiments. 
It is well known that those systems can exhibit an instability towards 
a regime where the cantilever settles into self-sustained oscillations. In 
this paper, we briefly review the classical theory of the optomechanical 
instability, and then discuss the features arising in the quantum regime. 
We solve numerically a full quantum master equation for the coupled sys- 
tem, and use it to analyze the photon number, the cantilever's mechanical 
energy, the phonon probability distribution and the mechanical Wigner 
density, as a function of experimentally accessible control parameters. We 
observe and discuss the quantum-to-classical transition as a function of a 
suitable dimensionless quantum parameter. 

1 Introduction 

Light interacting with matter can not only be scattered, absorbed and emitted 
by individual atoms, but it can also lead to mechanical effects. The radiation 
pressure of light was first directly observed in the seminal experiments of Nichols 
and Hull in 1901 and, independently, by Lebedev, where it exerted a torque on 
a pair of glass mirrors inside an evacuated chamber. Radiation pressure can also 
deflect the tail of comets (as first hypothesized by Johannes Kepler) , or change 
the path of asteroids. The mechanical effects of light become most pronounced 
in an optical cavity where the light intensity is resonantly enhanced, and where 
one of the end-mirrors is made movable, e.g., by being attached to a cantilever 
(Fig. [I). The pioneering theoretical and experimental works in this domain are 
due to Braginsky [TJ [2], In the 80s, strong effects were observed by the MPQ 
group of H. Walthcr in a setup using a macroscopic mirror [3]. More recently, 
the trend has been to exploit the tools of microfabrication to fabricate small 
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Figure 1: The basic optomechanical setup. 



cantilevers, nanobeams or other mechanical elements that can be affected by 
light. The small masses, high mechanical quality factors and (in some of the 
experiments) high optical finesse in these setups increase the optomechanical 
effects. Several recent experiments have demonstrated optomechanical cooling 
[41 El El O El HQl [11] , where the time-delayed light-induced forces lead to 
additional damping. Such a scheme may ultimately be employed to cool down 
to the ground-state of mechanical motion [TU [T3] ■ 

On the other hand, the light-induced forces can also lead to a negative con- 
tribution to the overall damping rate. At first, this increases the mechanical Q 
of the mechanical degree of freedom (hereafter simply referred to as the can- 
tilever), and can thus serve to amplify the response to any noise source acting 
on the cantilever. Once the overall "damping rate" becomes negative, which can 
happen simply by increasing the light power entering the cavity, we no longer 
have damping but instead an instability O [T3 EH E] . The cantilever starts to 
oscillate at its cigenfrequcncy, with an amplitude that at first increases exponen- 
tially and then saturates to a steady-state value. These self-induced oscillations 
have by now been observed in several experiments [TSJ Q33 HH1 121] ■ Theoretical 
studies predict an intricate attractor diagram |17j . which can display multiple 
stable attractors (i.e. possible oscillation amplitudes) for a given set of fixed 
external parameters. This attractor diagram has recently been observed and 
studied systematically in a low-finesse setup dominated by bolometric forces, 
which exhibited the unexpected feature of simultaneous excitation of several 
mechanical modes [2"T] . 

Similar physics has by now been observed in a variety of other systems which 
do not contain any optical elements. This includes driven LC circuits coupled 
to cantilevers [22], single-electron transistors and microwave cavities coupled to 
nanobeams [131 [Ml HS1 HS1 E3 HH [23 EO] , as well as clouds of cold atoms in an 
optical lattice inside a cavity [3TJ [35] . 

The common characteristic of all of these systems is that they contain some 
driven resonant quantum system (optical or microwave cavity, LC circuit, su- 
perconducting single-electron transistor), whose resonance frequency depends 
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on the motion of a mechanical degree of freedom (cantilever, nanobeam, defor- 
mation of a microtoroidal optical resonator, collective coordinate of a cloud of 
atoms) . Their Hamiltonian thus is typically of the form 

H = hujfi(x) a* a + hu>M &c + . . . , (1) 

where a is the photon annihilation operator for the driven resonator, whose fre- 
quency depends on the coordinate x = xzpf(c+ tf) of the mechanical oscillator 
with frequency u>m- Usually, it is possible to a very good approximation to 
expand uji to linear order in x, which is the case we will assume in the follow- 
ing. The additional terms not displayed in Eq. ([!} then describe the driving, as 
well as the damping and the fluctuation terms coupling to both oscillators (see 
below) . 

Cooling to the mechanical ground-state will generate the opportunity to 
observe a variety of quantum effects in such systems, including "cat" states [55] . 
entanglement [3H ISi] an d Fock state detection [TTJ [3j5] . For a recent review of 
optomechanical systems, see [37] > an d [3S] for the quantum noise approach to 
cooling. 

One question that may be asked about the quantum regime of these devices, 
that is now being approached experimentally, is how the instability discussed 
above changes due to quantum effects. We will answer this question partially 
in the present paper. We note that there have recently been some discussions 
of the quantum dynamics for the related instability in electronic systems in the 
literature [27j [28 29, 39 . The most important dimensionless parameter entering 
our analysis will be the "quantum parameter" 

C-^^, (2) 

a^FWHM 

which denotes the ratio between the mechanical zero-point fluctuation ampli- 
tude (a quantum parameter, tx Vh) and the width of the optical resonance (a 
classical quantity), measured in terms of displacement. After some rearrange- 
ment, one can see that this is essentially the "granularity" parameter employed 
in the discussion of the cold-atom experiment |32j . There, it was introduced 
by considering the total momentum kick a single photon would impart to the 
mechanical clement, as it is reflected multiple times before leaving the cavity. 
The granularity parameter then would be derived from the ratio of this kick to 
the mechanical momentum ground-state uncertainty. 

Increasing the quantum parameter £ will enhance quantum effects on the 
motion of the cantilever. These include the effects of the photon shot noise, 
as well as the mechanical zero-point fluctuations. The purpose of the present 
paper is to discuss these features in their dependence on £. We note that £ is 
rather small in the current optomechanical experiments (reaching up to about 
£ ~ 1CP 3 in |9J). However, given the large variety of analogous systems that are 
now being considered, we feel that it is justified to illustrate some of the salient 
features of the "quantum-to-classical crossover" by also analysing the regime 
£-1. 
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The remainder of this paper is organized as follows: We will first introduce 
the model Hamiltonian, and, in particular, discuss the full set of dimensionlcss 
parameters that are needed in our analysis. We then review the classical regime 
of self-induced oscillations. In particular, we will discuss the attractor diagram 
for the "resolved sideband regime" lom 3> k, which has not been discussed before 
in this context. This regime is currently of considerable interest, as it is crucial 
for ground-state cooling [T2 HH], and ujm/k ~ 20 has recently been realized 
experimentally [40 . Then we will turn to the full quantum model, that is first 
being discussed in terms of the rate equations that can yield the behaviour 
below the instability threshold. Afterwards, the full-blown quantum dynamics 
of self-induced oscillations will be analyzed using a master equation approach 
applied to the coupled system consisting of optical mode and cantilever. We will 
illustrate how the average mechanical energy, as a function of laser detuning, 
approaches the known classical result when the quantum parameter is sent to 
zero. In our present analysis, we focus on the steady-state and also assume a 
thermal bath temperature of T = 0. In real optomechanical experiments, one 
would presumably first cool down using the light field and then observe the 
onset of nonlinear dynamics as the detuning is varied. Investigation of these 
effects would entail studying the complete nonequilibrium time-evolution. 

2 Model and parameters 

In this section we present the Hamiltonian of the coupled cavity-cantilever sys- 
tem. A reduced set of dimensionlcss parameters determining the dynamics of the 
coupled system is identified. In particular, we introduce a quantum parameter, 
which is absent in descriptions of the classical dynamics (T7J [ST] and governs the 
crossover from classical to quantum behavior of the coupled dynamics of cavity 
and cantilever. 

2.1 Hamiltonian 

To describe a system of a mechanical cantilever coupled to a driven cavity, we 
consider the Hamiltonian 

H = H(-A - g(c+ c f )) a f a + hui M ctc + ha L (a + a f ) + H K + H r , (3) 

which is written in the rotating frame of the driving laser field of frequency 
lol, with an amplitude set by a^. The laser is detuned by A = u>l — w cav 
with respect to the optical cavity mode, described by photon annihilation and 
creation operators a and , and a photon number n cav = a' a. The cantilever 
(or, in general, mechanical element) of frequency % and mass m has a phonon 
number tim = c^c, and its displacement is given as x — xzpf(c + &), with the 
ground state position uncertainty (mechanical zero-point fluctuations) Xzpf = 
yjh/ (2tocjm)- The optomechanical coupling, between the optical field and the 
mechanical displacement, is characterized by the parameter g. In the simplest 
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case, with a movable, fully reflecting mirror at one end of an optical cavity, we 
have g = uj cav xzpF/ L, and thus g(c + & ) = oj cav x/L, with a radiation pressure 
force equal to i^ ra< j = a)ahg/xzPF- The decay of a photon and the mechanical 
damping of the cantilever are captured by H K and Hp M , respectively. They 
describe coupling to a bath leading to a cavity damping rate k and mechanical 
damping Tm- Note that each of the parameters A, g, ujm, ol has the dimension 
of a frequency 

2.2 Reduction to a set of dimensionless and independent 
parameters 

We now identify the dimensionless parameters the system dynamics depends on. 
Expressed in terms of the mechanical oscillator frequency ujm, the parameters 
describing the classical system are 



mechanical clamping 
cavity decay 
detuning 
driving strength 



k/um 
A/lj m 

V = 8\o L \ 2 g 2 /u M = oj cm n 2 E^J^ M mL 2 ). 



Here ££fax i s the light energy circulating inside the cavity when the laser is in 
resonance with the optical mode. The quantum mechanical nature of the system 
is described by the "quantum parameter" £, comparing the magnitude of the 
cantilever's zero-point fluctuations, xzpf, with the full width at half maximum 
(FWHM) of the cavity (translated into a cantilever displacement xfwhm) 

, xzpf 9 
quantum parameter : Q = = — . 

^FWHM « 

The resonance width of the cavity can be expressed as xfw hm — KL/w cav , 
where L is the cavity's length. The quantum parameter £ vanishes in the clas- 
sical limit h — > 0, as the zero-point fluctuations xzpf of the cantilever go to 
zero. The magnitude of £ determines the effect of quantum fluctuations on the 
dynamics of the coupled cavity-cantilever system. 

In later parts of this paper we will discuss the motion of the mechanical 
cantilever due to the driving of the cavity, both in a classical and a quantum 
mechanical picture. Such motion can be characterized by its energy E M , which 
in the classical case directly follows from the oscillation amplitude A of the can- 
tilever: Em, ci = \ muj M^ 2 ■ I n a quantum mechanical treatment, the energy is 
obtained from the expectation value of the occupation number of the oscillator: 
Em, qm = fi^M(^M), where we exclude the zero-point energy The dimensionless 
ratio of the cantilever energy Em to a characteristic classical energy scale of the 
system is then easily compared for the two approaches. To set this characteristic 
energy scale, we take the energy Eq = I^i^m^fwhm associated with an oscilla- 
tion amplitude xfwhm of the mechanical cantilever which moves the cavity just 
out of its resonance. Note that Em/E = (^4/xfwhm) 2 in the classical case, 
and Em/Eq = AC, 2 {fiu) in the quantum version. 
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3 Dynamics of the system 



In this section, we will first briefly recapitulate the results of a classical treat- 
ment of the optomechanical system |17[ 12 1] . We will concentrate on the regime 
of blue-detuned excitation of the cavity, where the cantilever motion is am- 
plified and self-induced oscillations can occur. Amplification behaviour of the 
coupled system (away from the regime of self-induced oscillations) can be un- 
derstood within a simple rate equation approach, which captures the effect of 
photon shot noise leading to fluctuations of the radiation pressure force acting 
on the cantilever |12j . The full quantum mechanical treatment, employing the 
numerical solution of a quantum master equation, can describe the crossover 
from heating/amplification to classical self-induced oscillations of the coupled 
system, where the quantum parameter £ = Xzpf/^fwhm governs the quantum- 
to-classical transition. 



3.1 Classical solution 

Heisenberg equations of motion for the cavity operator a and the cantilever 
position operator x can easily be derived from the Hamiltonian, Eq. [3] To 
investigate the purely classical dynamics of the coupled cavity-cantilever system, 
we replace the operator a(t) by the complex light amplitude a(t) and the position 
operator of the cantilever x by its classical counterpart. We thus arrive at: 

a = [i(A + g ) - -]a-ia L (4) 

XZPF £ 



Hg 



J M- h T — 

mxzpF 



a 



2 



- T M i ■ (5) 



Here fluctuations (both the photon shot noise as well as intrinsic mechan- 
ical thermal fluctuations) have been neglected, to obtain the purely deter- 
ministic classical solution. The variables t, x and a can be rescaled |17j as 
t = LUnjt; a = iotLOM /{Zocl); x — gx/(u)M%zPF) > so that the coupled equa- 
tions of motion contain only the dimcnsionless parameters V, A/w^, k/u)m, 
and T M /u M : 



da . . A _ . Ik. 1 

— - = M ha;) ]a+- 

dt lum 2u>m 2 

d 2 x , „ |2 T A/ di 

— =- = —x + V \a\ -= . 

dt 2 uj m dt 

Crucially, the quantum parameter £ cannot and does not feature in these equa- 
tions. 

Apart from a static solution x[t) = const, this system of coupled differential 
equations can show self-induced oscillations. In such solutions, the cantilever 
conducts an approximately sinusoidal oscillation at its unperturbed frequency, 
x(t) w x + Acos(tOMt). The radiation pressure affects the cantilever motion 
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rather weakly, so that the oscillation amplitude A varies only slowly and can 
be taken as constant during one oscillation period. The light amplitude then 
shows the dynamics of a damped, driven oscillator, which is swept through its 
resonance, see Eq. Q; an exact solution for the light amplitude a(t) can be 
given as a Fourier series containing harmonics of the cantilever frequency u>m 
HZ]: 



Kt)\ = 



with 



(6) 



1 



Jn(-A) 



2 in 



• (7) 

k/(2u m ) - i(x + A/ u>m ) 

The dependence of oscillation amplitude, A, and average cantilever posi- 
tion, x, on the dimensionless system parameters can be found by two balance 
conditions: Firstly, the total force on the cantilever has to vanish on average, 
and, secondly, the power input into the mechanical oscillator by the radiation 
pressure on average has to equal the friction loss. 

The force balance condition determines the average position of the oscillator, 
yielding an implicit equation for x, 



(x)=0 



J M-< 



(Prad) — 



mx Z pF 



(8) 



where the average radiation force, (P ra d) is a function of the parameters x and 
A. 

The balance between the mechanical power gain due to the light-induced 
force, P ra d = (-Fi-adi), and the frictional loss Pf r i c = Tm (x 2 ) follows from 



(xx)=0 {F xad x) = T M (x 2 }. 



(9) 



For each value of the oscillation amplitude A we can now plot the ratio between 
radiation power input and friction loss, P ra d/Pfric — (Pradi)/(rM(i 2 )), after 
eliminating x, using Eq. [8] This is shown in Fig. [2] Power balance is fulfilled if 
this ratio is one, corresponding to the contour line P ra d /Pfric = 1- If the power 
input into the cantilever by radiation pressure is larger than frictional losses 
(i.e., for a ratio larger than one), the amplitude of oscillations will increase, 
otherwise it will decrease. Stable solutions (dynamical attractors) are therefore 
given by that part of the contour line where the ratio decreases with increasing 
oscillation amplitude (energy) , as shown in Fig. [2] 

Changing the (dimensionless) mechanical damping rate Tm/^m will scale 
the plot in Fig.|2]along the vertical axis, so that the horizontal cut at one yields a 
different contour line of stable solutions [a changed input power V gives a similar 
scaling, but leads to further changes in the solution, as V also enters the force 
balance condition, Eq. (181]. Decreasing mechanical damping or increasing the 
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power input will increase the plot height in Fig.[2j so that the amplitude/energy 
of oscillation of the stable solution increases. 

While the surface or contour plots in Fig. [2] allow a discussion of general 
features of the self-induced oscillations, such as the multistabilities discussed 
in Ref. [T7], a slightly different representation of the classical solution is more 
amenable to an easier understanding of the particular dynamics of the system 
for a certain set of fixed system parameters. Figure [3]shows the cantilever energy 
Em,c\ — \ mL0 \i J ^ 2 m terms of the classical energy scale Eg = |ww^ f x| WHM as 
function of driving V and detuning A/lum- These are the parameters that can 
typically be varied in a given experimental setup. 

For sufficiently strong driving, self-induced oscillations appear around integer 
multiples of the cantilever frequency, A ps tilum- For a cavity decay rate k = 
0.5ojm assumed in Fig.[3j the different bands are distinguishable at lower driving; 
for larger k (or for stronger driving), the various 'sidebands' merge. For the 
lower-order sidebands, the nonzero amplitude solution connects continuously to 
the zero amplitude solution, which becomes unstable. This is an example of a 
(supercritical) Hopf bifurcation into a limit cycle. 

The vertical faces, shown gray in Fig. [3j for A sa 2u>m and A w 3ojm are 
connected to the sudden appearance of attractors with a finite amplitude. For 
example, while approaching the detuning of A = Ilom at fixed V (the solid 
line in Fig. prefers to V = 1.47 • 10~ 3 ), a finite amplitude solution appears, 
although A — remains stable. In Ref. [T7] the existence of higher- amplitude 
stable attractors and, correspondingly, dynamic multistability were discussed. 
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Figure 2: Classical self-induced oscillations of the coupled cavity-cantilever sys- 
tem. The radiation pressure acting on the cantilever provides an average me- 
chanical power input of P ra d. The ratio P-ad/Pfric of this power P ra( j vs. the loss 
due to mechanical friction, Pf r i C , is shown as a function of the detuning A and the 
cantilever's oscillation energy Em, at fixed laser input power V . The oscillation 
energy Em = mu! M A 2 /2 is shown in units of -Bo, where Em/E = (A/xfwhm) 2 - 
Self-induced oscillations require P ra d = Pmc- This condition is fulfilled along 
the horizontal cut at P ra d/Pric = 1 (see black line and the inset depicting the 
same plot, viewed from above). These solutions are stable if the ratio P-ad/Pric 
decreases with increasing oscillation amplitude A. The blue regions at the floor 
of the plot indicate that P ra( j is negative, resulting in cooling. The cavity decay 
rate is k = 0.5u>m, the mechanical damping is chosen as Ym/^m = 1-47 • 1CP 3 , 
and the input power as V = 6.05 • 1CP 3 ; these parameters are also used in 
figures [3] [4j [5j and[6j and will be referred to as r^ f and V* . 
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cantilever energy 




Figure 3: Cantilever oscillation energy E M oc A 2 versus detuning A and laser 
input power V. This plot (in contrast to Fig. [2]) shows only the stable oscilation 
amplitude, but as a function of variable input power. The particular value 
V* corresponding to Fig. [2j and the resulting profile of oscillation amplitudes 
are indicated by a black line. The green floor of the plot indicates regions 
without self-induced oscillations. The other system parameters are as in Fig. [2] 
The continuous onset of the self-oscillations in the sidebands at A/u>m = 0, 1 
(which merge for the present parameter values) represents a supercritical Hopf 
bifurcation, from A — to A ^ 0. At higher sidebands, an attractor with a 
finite A ^ appears discontinuously, while A = remains a stable solution. 

3.2 Rate equation approach 

Before embarking on a full quantum- mechanical treatment of the coupled cavity- 
cantilever system, it is instructive to discuss a more simple method to capture 
some nonclassical effects, in particular the response of the cantilever to the 
photon shot noise. For that purpose, we consider the shot noise spectrum of the 
driven cavity, decoupled from the cantilever, 

S ff (uj) = ) S nn (uj) = ) n- — , (10) 

\x ZPF J \xzpfJ (^ + A) 2 + {k/2) 2 
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where 

8( 2 (A/ujm) 2 + (k/2ujm) 2 { 1 

is the mean number of photons in the cavity. The maximum occupation n max = 
Voj m /(2k 4 £ 2 ) = Aa 2 L / k 2 occurs at zero detuning. We note that in using the 
unperturbed, intrinsic shot noise spectrum for an optical cavity in the absence 
of optomechanical effects, we neglect the modification of that spectrum due to 
the backaction of the cantilever motion. 

The asymmetry of the shot noise spectrum is important for the dynamics of 
the cantilever. The spectral density of the radiation-pressure force at positive 
frequency lum (negative frequency — u>m) yields the probability of the cavity 
absorbing a phonon from (emitting a phonon into) the cantilever |12j . 

For a red-detuned laser impinging on the cavity (A < 0), the cavity's noise 
spectrum peaks at positive frequencies and the cavity tends to rather absorb 
energy from the cantilever. As a consequence, the mechanical damping rate for 
the cantilever is increased, leading to cooling if one starts with a sufficiently 
hot cantilever. In the opposite Raman-like process taking place at A > 0, 
a blue-detuned laser beam will preferentially lose energy to the cantilever, so 
that it matches the cavity's resonance frequency. The effective optomechanical 
damping rate, 

r op t = C 2 K 2 [S nn (+ujM) - S„„(-wm)] , (12) 

is then negative. The corresponding heating of the mechanical cantilever is 
counteracted by the mechanical damping Tm ■ Simple rate equations for the 
occupancy of the cantilever yield a thermal distribution for the cantilever phonon 
occupation number um, with |12j 

(c T c) = (n M ) = — . (13) 

1 opt + 1 M 

The effective temperature, T e g, is related by (hj^+1) /(tim) — exp[tkJM /(fes^eff)] 
to the mean occupation number. The equilibrium mechanical mode occupation 
number, n t h, is determined by the mechanical bath temperature, which is taken 
as zero in the following. In contrast to first appearance, the mean occupation 



number of the cantilever given in Eq. (13 1 does not depend on the quantum 



parameter £, as C Snn is independent of £. This is because S nn ~ n ~ 1/Q 



2 



see Eq. (111. The cantilever energy, therefore, only trivially depends on the 
quantum parameter as Em/Eq = 4£ 2 (hm), so that it vanishes in the classical 
limit, where C 2 cx /i — > 0. 



In general, the phonon number in Eq. ( 13 1 can increase due to two distinct 
physical effects: On the one hand, the numerator can become larger, due to the 
influence of photon shot noise impinging on the cantilever, represented by S nn . 
On the other hand, the denominator can become smaller due to r opt becoming 
negative. In the latter case, the fluctuations acting on the cantilever (both 
thermal and shot noise) are amplified. This effect is particularly pronounced 
just below the threshold of instability, where Tm + Topt = (see below). 
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In the resolved sideband limit k « (at weak driving) the cantilever 
occupation (um) will peak around zero detuning, where the number of photons 
in the cavity is large, and around a detuning of A = u>m- At the latter value of 
detuning the aforementioned Raman process is maximally efficient as a photon 
entering the cavity will exactly match the resonance frequency after exciting 
a phonon in the cantilever. This dependence of cantilever occupation (or the 
corresponding energy) on the detuning is shown in Fig. [4] 

The approach sketched above can be modified slightly to take account of the 
modification of the cavity length due to a static shift of the cantilever mirror 
by radiation pressure. Approaching the resonance of the cavity from below, the 
increasing number of photons inside the cavity will increase the cavity length due 
to their radiation pressure on the mirror, bringing the system even closer to the 
resonance. The effect of the static shift of the mirror on the mean occupation 
of the cavity can be included self-consistcntly, leading to the tilt of the peak 
around the resonance, shown by the dash-dotted line in Fig. |4|a). The same 
figure also includes results of the full quantum mechanical approach, which will 
be discussed in the next section. 

For larger n, the two peaks in the cantilever excitation merge. Higher-order 
sidebands are not resolved within this approach, since they would require taking 
care of the modification of Sff due to the cantilever's motion. 

Classical self-induced oscillations occur in a regime of larger driving, where 
the optomechanical damping rate r opt of Eq. (12) becomes negative. They 
appear once amplification exceeds intrinsic damping, i.e. when r opt + Tm < 0. 
The simple rate equation approach lacks any feedback mechanism to stop the 
divergence of the phonon number. The classical solution demonstrates how this 
feedback (i.e. the resulting change in the dynamics of the radiation field) makes 
the mechanical oscillation amplitude saturate at a finite level. In addition, it 
shows the onset of self-induced oscillations to occur at a smaller detuning, due 
to the effective shift of the cantilever position explained above. 
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Figure 4: Cantilever energy versus detuning for a cavity driven below [(a),(c)] 
and above [(b), (d)] the onset of self-induced oscillations. Note Em/Eq = 
AQ 2 {um)- (a) Below the onset, the cantilever amplitude would vanish ac- 
cording to the classical analysis that does not incorporate fluctuations. How- 
ever, the cantilever is susceptible to the photon shot noise (the parameters are 
k/lu m = 0.1, V = 8.4 ■ 1CT 3 , T m /lu m = 5 • 1(T 3 , and C = 1.0), leading to finite 
phonon numbers in the cantilever, particularly around the resonance A = 
and at the first sideband A = lom (see main text). This is captured by the 
full quantum master equation, as well as (approximately) by the rate equation, 
whose results improve when taking into account the corrections due to the shift 
of the cantilever position x. (b) For stronger driving, the classical solution yields 
self-oscillations (the parameters are V* , T* M as in Fig. [3] but k/l>m = 0.3). The 
rate equation correctly predicts the onset of the linear instability, but not the 
nonlinear regime. [The shift in x was not taken into account, hence the slight 
discrepancy vs. the classical solution] The master equation results are shifted 
to lower detuning and describe sub-threshold amplification and heating as well 
as self-induced oscillations above threshold, modified and smeared due to quan- 
tum effects (as shown for a quantum parameter of £ = Jzpf/^fwhm = 1)- (c) 
Including the zero-point fluctuations in a semiclassical approach via Langevin 
equations gives results that agree well with both the results from the rate equa- 
tion and the full master equation, shown here for parameters as in (a), (d) 
Above the onset of self-induced oscillations the semiclassical approach mimicks 
results from the quantum master equation partially. The parameters for this 
plot are k/lu m = 0.3, T M = 50T* M , V = 20V*, ( = 1. 



13 



In Fig.[i|b) we show results for the detuning dependence of the mean energy 
of the cantilever obtained from this rate equation approach below the thresh- 
old of classical self-induced oscillations. The coupled cavity-cantilever system 
acts as an amplifier of fluctuations, increasing the occupation of higher number 
states of the cantilever well before classical oscillations set in. At the onset 
of classical self-induced oscillations the rate equation result diverges. A full 
quantum-mechanical treatment describes the crossover of the cantilever dynam- 
ics from quantum-fluctuation induced heating to self-induced oscillations as will 
be discussed now. 



3.3 Quantum master equation method 

The evolution of the coupled quantum system comprised of the cantilever and 
the optical cavity is described by the Hamiltonian of Eq. ^ . Dissipation arises 
from the coupling of the mechanical mode to a bath and due to the opening 
of the cavity to the outside. While the former results in mechanical damping 
with a rate Tm, the latter is associated with the ringdown rate of the cavity 
K. In the present paper, we will assume the mechanical bath to be at zero 
temperature, where quantum effects are most pronounced in steady state. A 
future, more realistic treatment, should relax this assumption and deal with the 
noncquilibrium dynamics that results when a mechanical system is first cooled 
optomechanically and then switched to the unstable side. 

The system can be described by a reduced density matrix p for the mechan- 
ical cantilever mode and the optical mode of the cavity. In the frame rotating 
at the laser frequency, the time evolution of the density matrix p is given by 

at in 

where T>[A] = ApA^ — ^A^Ap — ^pA^A denotes the standard Lindblad opera- 
tor. The Hamilton operator Hq describes the coherent part of the evolution of 
the coupled cavity-cantilever system, 

H = Hn + H K + H v . 



For the numerical evaluation, we rewrite Eq. 14 as dp/dt — Cp, with a Liouvillian 
superoperator C. We then interpret the density matrix as a vector, whose time 
evolution is governed by the matrix C The density matrix at long times (in 
steady state) is then given by the eigenvector of C with eigenvalue 0. The 
numerical calculation of this eigenvector is much more efficient than a simulation 
of the full time evolution. Since we are dealing with large sparse matrices, it 
is convenient to employ an Arnoldi method that finds a few eigenvalues and 
eigenvectors of C by iterative projection. For Hermitean matrices, the Arnoldi 
method is also known as the Lanczos algorithm. 

In practice, the numerical approach used here sets strong limits on the di- 
mension of the Hilbert space. We need to take into account the N a lowest Fock 
states of the cavity and the N c lowest Fock states of the mechanical cantilever, 
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resulting in a Liouvillian superoperator with (N a x N c ) 4 elements. This puts 
more severe restrictions on our treatment of the coupled cavity-cantilver system 
than encountered in similar treatments of comparable systems. For example, 
nanoelectromechanical systems, where an oscillator is coupled to a normal-state 
or superconducting single-electron transistor (SET), will have to account for 
only a very limited number of charge states of the SET (namely those few in- 
volved in the relevant transport cycle). As a consequence, a larger number of 
Fock states can be included, e.g., 70 number states of the oscillator were kept in 
Ref. [28]. In some cases it was furthermore considered sufficient to treat only 
the incoherent dynamics of the mechanical oscillator, i.e., only the elements of 
the density matrix diagonal in the oscillator's Fock space, thereby reaching 200 
number states of a mechanical mode coupled to a normal-state SET [IT]. The 
restricted number of Fock states that can be considered here makes it more dif- 
ficult to fully bridge the gulf to the classical regime of motion of the mechanical 
cantilever. [(N a , N c ) = (8,16) for Fig. ga),(c),(d), (4,22) for Figs. Qb), [5] and 
for the first two panels of [6] (3,35) for the last panel of Fig. [6] 

A first comparison of results of the quantum master equation to the classical 
solution and the results of the rate equation was already shown in Fig. [4] We find 
that the full quantum results do not qualitatively differ from the rate equation 
results provided the parameters are chosen sufficiently far from the onset of self- 
induced oscillations. However, only the quantum master equation approach is 
able to describe the crossover from sub-threshold quantum fluctuations (where 
Em oc £ 2 ) to the large classical cantilever energies associated with self-induced 
oscillations. 

In Fig. [5] we demonstrate the influence of the quantum parameter ( = 
^zpf/^fwhm governing the crossover between the classical and the quantum 
regime. 

Figure |5ja) shows the cavity photon number, normalized to its value at res- 
onance, n max . For our choice of driving parameter V, the maximal occupation 
n max is low, so that a small number of Fock states suffices for describing the cav- 
ity in the quantum master equation. This allows to account for enough number 
states of the cantilever to reach the regime of self-induced oscillations. The clas- 
sical solution (solid black line) consists of the broad Lorentzian of the isolated 
cavity, on top of which additional peaks appear. These are due to the classical 
self-induced oscillations occuring at the sidebands A = ujm > 2u>m , ... in the cou- 
pled cavity-cantilever system. Figure[5jc) displays the cantilever energy E M /E 
as a function of the detuning, A/ujm, with features that parallel those found 
for the photon number. The classical curve in (b), shown in black, corresponds 
to the cut indicated by the solid line in Fig. [3] For the chosen driving power, 
the second sideband at A = u>m just starts to appear, while the first sideband 
is merged with the resonance at A = 0, which shows up as a slight shoulder. 
The sharpness and strength of these features also depend on the values of me- 
chanical damping and cavity decay rate. Results of our solution of the quantum 
master equation are shown for three different values of the quantum parameter 
C = ^zpf/^fwhm- Due to restrictions of the numerical resources, it was not 
feasible to map out a wider range of values of the parameter £, although the 
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range analysed here already suffices to describe the quantum-classical crossover. 

The quantum master equation shows results that are qualitatively similar to 
the classical solution in the regime of self-induced oscillations, with the peaks 
being progressively broadened, reduced in height, and shifted to lower detuning 
for increasing values of the quantum parameter £. Numerical evidence indicates 
that quantum correlations between the cantilever position operator x and the 
photon operators a,', a may cause the observed shift. As expected, the dis- 
crepancy between the quantum mechanical and the classical result reduces with 
diminishing quantum parameter £. In Fig. [5jb) , we show the dependence of 
the cantilever energy on the quantum parameter, for two different values of the 
detuning. In the sub-threshold regime of amplification/heating the cantilever 
energy scales as £ 2 , as discussed above. In any case, the classical limit is clearly 
reached as £ — > 0. 

At the second sideband a classical solution of finite amplitude coexists with 
a stable zero-amplitude solution (compare Fig. [2] and last panel of Fig. [6]). The 
black curve in Fig. [5jb) , showing the finite amplitude solution, may therefore 
deviate substantially from the h — > limit of the quantum mechanical result. 
In general, the average value of Em, shown here, will be determined by the 
relative weight of the two solutions (which are connected by tunneling due to 
fluctuations), as well as fluctuations of Em for each of those two attractors. 

3.4 Langevin equation 

To get an estimate of the influence of quantum fluctuations, we compare the 
results of the quantum master equation to numerical simulations of classical 
Langevin equations that try to mimick the quantum noise. The resulting de- 
scription of the quantum-to-semiclassical crossover is illustrated in Figs. Qc). 
To imitate both the zero-point fluctuations of the mechanical oscillator and the 
shot-noise inside the cavity, we add white noise terms to Eqs|4jind|5] 

d = [i(A + g ) - -] a - iai + \J k/2 a in (15) 
X Z PF 2. 

x=—cj m x-\ — \a\ 2 — T M x + JfiwMr/m^, (16) 

mx Z pF 

where (a m ) = (£) = and (a in = = *(* ~ 0- The co- 

efficients in front of the noise terms are chosen such that in the absence of 
optomechanical coupling we obtain the zero-point fluctuations, i.e. (|a| 2 ) = 0.5 

away from resonance and mu ^ M (x 2 ) — The mean zero-point energy of the 

cantilever is substracted from the curve in FigQc). 

For parameters below the onset of self-sustained oscillations, this semiclassi- 
cal approach leads to good qualitative agreement with the quantum mechanical 
description, as can be seen in Fig. |4jc) for parameters that are the same as 
those of |4^a). Still, the Langevin approach can mimick the results from the 
master equation only partially. In particular, the approximation gets worse 
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when dealing with low photon numbers. This is because the Langevin equation 
introduces artificial fluctuations of the radiation pressure force in the vacuum 
state. Indeed, |a| 2 has a finite variance even in the ground state of the photon 
field, in contrast to a) a. 



3.5 Wigner density and phonon number distribution 

In figure [6] we go beyond the average cantilever phonon number and present 
results both for the phonon number probability distribution, as well as the full 
Wigner density of the cantilever, defined as 

1 f + °° 

W(x,p) = — (x-y\p\x + y) e 2l ™' h dy. (17) 

7™ J-oo 

This figure demonstrates the different nature of the cantilever dynamics in 
the sub-threshold regime and above threshold, where self-induced oscillations oc- 
cur. Below the threshold (for a detuning A a = — 0.45wm as indicated in Fig. [5] 
quantum parameter £ = 1, and other parameters as in Fig. [5]) the occupation of 
the cantilever is thermal, with an effective temperature determined by the effec- 



tive optomechanical and mechanical damping rates, cf. Eq. (13). Consequently, 
the Wigner density shows a broad peak around the origin of the x — p plane of 
cantilever position and momentum (the static shift of the cantilever due to the 
radiation pressure is very small). For a detuning of Af, = — 0.2um, self-induced 
oscillations occur. The probability distribution for the phonon number shows 
some thermal broadening, but an additional peak appears at a finite phonon 
number. In the Wigner density plot this results in a crater-like feature, which 
corresponds to a mixture of coherent states with essentially fixed amplitude 
but arbitrary phases. This captures the fact that the phase of the self-induced 
oscillations is completely arbitrary also in the classical solution. The energy 
corresponding to the phonon number at which the distribution peaks, compares 
fairly well to the oscillation energy obtained from the classical solution. Only 
the shift towards lower values of detuning as shown in Fig. (5^b) puts restrictions 
on a detailed quantitative comparison. 
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Figure 5: Comparison of classical and quantum results, (a) Number of photons 
inside the cavity as a function of detuning, and (c) energy of the cantilever versus 
detuning for Pj^, V* and k/um = 0.5. The dotted curves show results from 
the quantum master equation for different values of the quantum parameter 
C = 1.3 (pink) , £ = 1.0 (green) and £ = 0.7 (blue), which are compared 
with the solution of the classical equations of motion (black solid curve). As 
C — > 0, the qantum result approaches the classical curve. See main text for 
a detailed discussion, (b) The energy of the cantilever as a function of the 
quantum parameter £ for fixed detunings Aj/um = —0.2 and A c /uim = 0.4 
(the detuning value A a indicated in (b) is used in Fig. [fjj. (d) Fano factor 
((^m) — (^M) 2 )/(«Af) vs. detuning, for ( = 1. For a coherent state whose 
occupation number follows a Poisson distribution, the Fano factor is 1 (dashed 
black line). Close to the resonance (and far away from it, where (um) = 0), 
the results of the quantum master equation approach this value. The Fano 
factor becomes particularly large near the second sideband, where we observe 
coexistence of different oscillation amplitudes (see Fig. [6| . 
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Figure 6: Distribution functions P(jim) of the cantilever occupation and Wigner 
functions W{x,p) [rescaled by xzpfPzpf] of the cantilever for A a = — 0.45wm, 
At, = — 0.2wm, Ad = 1.72lom [corresponding to the detuning values also indi- 
cated in Fig. [5^b); further parameters as in Fig. [5] with C = 1-0; for A d the 
mechanical damping rate is reduced to Tm/<^m = 1-2 • 10 -3 ]. Below the thresh- 
old of self-induced oscillations, a broadened distribution is found corresponding 



to an increased effective temperature, cf. Eq. (13 1 (left panels, A Q ); self-induced 
oscillations are visible as a finite amplitude ring in the middle and the right 
panel. Dynamical multistability (i.e. co-existence of several attractors) in the 
classical solution becomes apparent both in the distribution and the Wigner 
density, where a double-peaked structure develops. 



For a value of the detuning located in the second sideband, A^ = 1.72u>m, we 
find a probability distribution with a peak for the occupation of the cantilever 
ground state, and a broader peak at a finite occupation number (mechanical 
damping is slightly decreased to display more pronounced features). Likewise, 
the Wigner density consists of a sharp peak at the origin, surrounded by a 
broader ring representing finite amplitude oscillations. This corresponds to 
the existence of two stable attractors in the classical analysis, with vanishing 
and finite oscillation amplitude, respectively. Similar results for the Wigner 
densities were found in Ref. [28 for a cantilever driven by a superconducting 
single-electron transistor. 



4 Conclusion 

We presented a fully quantum mechanical treatment of a driven optical cavity 
coupled to a mechanical cantilever by radiation pressure. Light-induced forces 
can yield a negative contribution to the damping of the cantilever, causing 
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amplification of fluctuations and even instabilities of the cantilever dynamics. 

In the present paper we first reviewed briefly the classical solution and dis- 
cussed the existence of self-induced oscillations and the resulting attractor di- 
agram of the system. We paid particular attention to the resolved-sideband 
regime k <C %, which is now increasingly studied in experimental setups. 
Here the instabilities clearly occur at sidebands, where the detuning matches 
an integer multiple of the mechanical frequency 

Within a simple rate equation approach, we were able to discuss the influ- 
ence of the photon shot noise and quantum fluctuations well below the instability 
threshold. The full quantum-mechanical treatment, based on a numerical solu- 
tion of the quantum master equation, is able to completely describe both regimes 
(below and above threshold) . It has been complemented by numerical studies of 
a Langevin equation that includes the zero-point fluctuations in a semiclassical 
way. We studied the crossover between the quantum and classical regime, which 
is governed by the quantum parameter, ( = ^zpf/^fwhm j denoting the ratio 
between the mechanical zero-point fluctuation amplitude and the width of the 
optical resonance. Signatures of the self-induced oscillations are also found in 
the full quantum mechanical solution, even at larger values of (. In regions of 
dynamical multistability, the different attractors show up simultaneously in the 
steady state of the cantilever, since the quantum noise can induce transitions 
between those attractors. Finally, we characterized the mechanical motion in 
the various regimes by discussing the phonon number probability distribution 
as well as the Wigner density. 
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